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The mass distribution of fragments that are ejected at a given velocity for impact 
craters is modeled to allow extrapolation of laboratory, field, and numerical results to 
large-scale planetary events. The model is semi-empirical in nature and is derived from 
(l) numerical calculations of cratering and the resultant mass versus ejection velocity, (2) 
observed ejecta blanket particle size distributions, (3) an empirical relationship between 
maximum ejecta fragment size and crater diameter, (4) measurements and theory of 
maximum ejecta size versus ejecta velocity, and (5) an assumption on the functional 
form for the distribution of fragments ejected at a given velocity. This model implies 
that for planetary impacts into competent rock, the distribution of fragments ejected at 
a given velocity is broad, e.g. 68% of the mass of the ejecta at a given velocity contains 
fragments having a mass less than 0.1 times a mass of the largest fragment moving at 
that velocity. The broad distribution suggests that in impact processes, additional com- 
minution of ejecta occurs after the upward initial shock has passed in the process of the 
ejecta velocity vector rotating from an initially downward orientation. This additional 
comminution produces the broader size distribution in impact ejecta as compared to that 
obtained in simple brittle failure experiments. 
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1.0 INTRODUCTION 
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Meteorites come from collisions on objects varying in size from asteroidal to terres- 
trial planets 1-8 . Massive extinctions on the earth may have resulted from large impacts 9 . 
The origin of meteorites is dependent on the maximum-size fragment that can be ejected 
at velocities exceeding the planetary escape velocity, whereas the biological consequence 
of large impacts is dependent on the number of mass of fine particles (diameters <10 -6 
m) that are ejected high into the atmosphere. 10 

A model was developed for the mass distributions of fragments that are ejected at a 
given velocity for impact events 10 . The model is semi-empirical in approach and is con- 
sistent with empirical data from small scale impact experiments, nuclear, and chemical 
explosion data, ejecta fragment size vs. velocity data, largest ejecta fragment size vs. 
crater size data and total mass ejected vs. impact energy. In addition, the model is 
based upon calculations of mass ejected as a function of ejection velocity and maximum 
fragment size as a function of strain rate models. Using the above data base, the distri- 
bution function for the fragments was derived. In this paper we have extended our pre- 
vious approach and have found additional data that strengthen those conclusions. 


2.0 CUMULATIVE MASS DISTRIBUTION OF FRAGMENTS 

The objective of this paper is to determine the distribution of fragments that are 
ejected at a given velocity during an impact event. Shown in figure 1 is a typical 
impact-induced flow field at a given time. The figure demonstrates the decrease in 
ejecta velocity, V, with decreasing radius, r, given approximately by V oc r . The crater 
plume consists of particles ejected at positions r at times t and the velocity of ejection 
varies as a function of position along the surface. Thus the plume necessarily changes 
shape during its evolution. 

By definition, the cumulative amount of mass M c of fragments ejected at all veloci- 
ties of mass greater than m is given by 

v 

Mc “ / f ( m ’ m bv( V )) dV (*) 

^itin 


where 


d M„ 
d V 


is the amount of mass ejected at velocities between V and V+dV, and 


f(m,mbv(V)) is the unknown distribution of fragments ejected at V. The expression for 
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the amount of mass ejected at a given velocity, M ev , can be determined from the fits to 
the cumulative amount of mass ejected at velocities greater than V obtained either from 
theoretical calculations 10-13 (see figure 2) or experimental data 14-15 . 


3M CV 
d V 


f Mt 




V min <LV <L V MAX 


( 2 ) 


where M T is the total amount of mass ejected, and V min and V max are the minimum and 
maximum velocities of ejection of solid particles. A discussion of the thermodynamic 
constraints on the maximum velocity of ejection is given by O’Keefe and Ahrens 10 . The 
exponent £ is related to the mass scaling parameter 16 , a, 

by f=6a/(a-3) (3) 

The mass scaling parameter specifies the relation between M c and m b , the maximum 
fragment mass. 

M 

= [l-(m/m b ) o/3 ] (4) 

The minimum velocity of ejection varies with the planetary strength at small scales and 
Rg at large scales, where R is the size of the impactor and g is the planetary gravity. 

The key assumption in the theory is that the functional form of the distribution of 
fragments ejected at a given velocity is similar to the distribution function of the frag- 
ments in the ejecta blanket. For explosive events Jaeger et al. 17 have found that this is 
the case. With this assumption the cumulative fraction of mass fragments of mass 
greater than m is given by 



where m bv is the mass of the largest fragment ejected at a given velocity. 

In addition, we assume that the mass of the largest fragment ejected at a given 
velocity is a function of the velocity of ejection and is given by 

= (V/V ml .)- ( (6) 

m b 

where 6 is an unknown parameter to be determined 
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The above functional form can be justified by comparison to explosion cratering 
data 18 , planetary cratering observational analysis 19 , and theoretical considerations. 
Ejecta velocities were measured on a 100 ton TNT above ground explosion called Middle 
Gust III, and are summarized by Schoutens 18 . The size-dependent ejection velocity limit 
inferred from film recording and analysis by Wisotski 20 and corrections for deceleration 
in air 21 is shown in figure 3. The range of launch velocities correlates with 6 ~ 3. 

The theory developed by Grady and Kipp 22 can be used to determine an expression 
for 6. They found that the mass of the largest fragment produced is a function of the 
strain-rate 

m b oc e 31/(1 +3) (7) 


where / ~ 8 for rocks. 

The strain-rate is related to the velocity field by 

av 


€ OC 


dr 


(8) 


where r is, again, the surface radial coordinate. The surface velocity field has the follow- 
ing form for impacts 16 . 


V oc 



( 9 ) 


where R is the radius of the impactor. Evaluating the surface strain-rate, the mass of 
the largest fragment created at a given velocity is 


m b oc 



fg+3)/ 

1+3 


( 10 ) 


and thus 

6 = (e+3)//(/+3). (ii) 

For a variety of materials, 6 ranges from 3 to 3.6. The expression given by equation 12 
for cumulative mass distribution of fragments in the ejecta blankets is compared to 
experimental data in figure 6. In the case shown, 6 = 3 and £ = 1.4. The only unk- 
nown parameter in equation 12 is /?, which is in the exponent in mass distribution of 
fragments ejected at a given velocity. Referring to figure 6, we get good agreement with 
the fragment size data when f) ~ 0.5. Note that this has the same value as the 
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cumulative mass distribution exponent a (Eq. 6). This means that the distribution of 
fragments ejected at a given velocity is polydispersed and that in the flow and ejection 
process, grinding and crushing has occurred since fracture processes are more mono- 
dispersed 22,26 . 


The integral for the cumulative mass distribution in the ejecta blanket can be 
evaluated by substituting equations 2, 5, and 6 in equation 1, which yields the following 



= 0 for m < m s 

where m s is the mass of the smallest fragments. Note that in the case of both impact 
craters 23 and explosion events 17 that there is a characteristic scale size in which the 
number of smaller fragments decreases. 

The above expression (equation 12) can be compared to experimental data on the 
cumulative mass distributions fragments ejected from impact craters 23-25 . Both impact 
and explosion fragments distributions can be described by expressions of the following 
form 11,26,27 . 

M c /M t = |l-(m/D, b )“ /3 l (13) 

where m b was found by Gault et al. 23 to correlate with the total mass ejected from the 
crater for large range of crater sizes. Shown in figure 5 is that data which is described 

by 

m b = 0.2Mj 8 (14) 

The following is an expression for the mass fraction of fragments that are ejected 
with velocities in fraction range AV at V. 


AM CV f 

AV 1 

i v ,-*«> 

1 __ 

( — ) 3 
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6a 

3 
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V * 
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( m b j 
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( 15 ) 
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3.0 CONCLUSIONS 

1. The ejecta velocity, V, of surface material achieved upon impact of a projectile 
of radius, R, at a radius, r, is approximately given by V cc (R/r) 3 . 

2. The cumulative mass of impact ejecta launched at a given velocity is described 
by a power law function, in which the exponent £ is related to the exponent, a, of the 
power law function which describes the cumulate mass, versus, mass greater or equal to 
a given value, of the total ejecta distribution by £=6c*/(a'-3). 

3. The velocity of the largest fragments in a given cratering ejecta distribution is 
described by a power law function of ejecta velocity where the scaling exponent, 6, is 
given by 6 = (£+3)/ /(/ +3) where l is the parameter in the expression for the largest 
fragment m b produced upon brittle failure at varying strain rates as described by Grady 
and Kipp 22 . The most massive fragment launched at a given velocity, m b , is found to 
depend on velocity as V, m b oc V -3 , in agreement with results of Vickery 19 . 

4. Both the cumulative mass of ejecta and the mass of ejecta launched to a velo- 
city, V, and the total mass greater than m in the ejecta have similar power law distribu- 
tion. 

5. Brittle fragments distributions obtained in simple failure experiments are very 
monodispersed, as to their size, in comparison to impact ejecta distributions which are 
very broadly dispersed with mass to sizes traveling at a given velocity. We infer that 
additional comminution occurs during impact cratering, presumably during the processes 
of ejecta rotation which take place on the particle velocity impacted by the shock and 
changes its direction from being primarily downward directed to achieving a net outward 
and upward motion. 
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Figure Captions 


Fig. 1. Cratering flow field at a time, t, given by d/Vt=2.0, where V is impact velo- 

city and d is projectile diameter for water projectile striking silicate half space 
at 5 km/sec. The relationship between the ejecta plume and surface velocity 
field is indicated. The velocity, V, of each fragment of ejecta is approximately 
related to its take-off radius, r, by V cc (R/r) 3 , where R is the projectile radius. 

Fig. 2. Cumulative ejecta mass at velocity less than V versus ejecta velocity. Calcu- 

lations for impacts of solid silicate and ice at 5 km/sec are indicated. Here 

M c /M t = M' c /M x -1. 

Fig. 3. Observed ejecta velocity, versus, ejecta diameter, A, 5 to 10 sec after explosion 
of MIDDLE GUST III, a 100 ton TNT surface explosion conducted over par- 
tially saturated Pierre shale 18 . Size dependent ejecta velocity limit correspond 
to 6 — 3.0. 

Fig. 4. Ejecta fragment size versus ejection velocity inferred from gravity scaling from 
Lunar (a) and (b), and Martian (c) craters 18 mass of fragments correlates with 
6 — 3 for relation m b oc V -3 shown. 

Fig. 5. Mass of the largest fragment versus total ejected mass for explosive and 
impact cratering events 23 . 

Fig. 6. Cumulative mass M c (ejected at velocity V) normalized by the total mass M T 
(ejected at velocity V) versus fragment size normalized to the largest fragment 
ejected at velocity V. Data is for a value of £ = 1.9 as indicated and theoreti- 
cal curves are present for 0 = 0.4, 0.5, and 0.6. Here M c /M x = M c '/M x ~l. 
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